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Abstract 

We have calculated electronic and optical properties of Si/BeSeo.4iTeo.59 heterostructures by 
a semiempirical sp 3 s* tight-binding method. Tight-binding parameters and band bowing of 
BeSeo.41Teo.59 are considered through a recent model for highly mismatched semiconductor al- 
loys. The band bowing and the measurements of conduction band offset lead to a type II het- 
erostucture for Si/BeSeo.4iTeo.59 with conduction band minimum in the Si layer and valence band 
maximum in the BeSeo.41Teo.59 layer. The electronic structure and optical properties of various 
(Si2)n/(BeSeo.4iTeo.59)-m [001] superlattices have been considered. Two bands of interface states 
were found within the bandgap of bulk Si. Our calculations indicate that the optical edges are 
below the fundamental bandgap of bulk Si and the transitions are optically allowed. 

PACS numbers: 73.21.Cd,73.21.Fg,78.67.De 



I. INTRODUCTION 



Silicon completely dominates the present semiconductor market with the greatest shares 
of all sales. It has an ideal bandgap (1.12 eV) for room temperature operation and its oxide 
(SiC>2) provides the necessary flexibility to fabricate millions of devices on a single chip. High 
integration implies high-speed operation that is limited by the interconnect propagation 
delay of the signal between devices. This constraint suggests that the integration of Si 
micro-electronics might be aided by optical interconnection. Unfortunately, silicon does not 
respond strongly to optical excitations because it is an indirect bandgap semiconductor: the 
band extrema for electrons and holes are located at different points in the Brillouin zone 
(Fig. Q). Therefore, intrinsic formation or recombination of electron-hole pairs becomes a 
three-particle event, which is weaker than a two-body process. Therefore, there are a number 
of areas where Si cannot compete and this has allowed other materials to dominate such as 
optoelectronics where the indirect bandgap precludes the use of Si to produce efficient light 
emitting diodes and lasers. If the performance of Si transistors or circuits could be improved 
by addition of another semiconductor material then numerous new applications could open 
up. Silicon-germanium is one such material which may be epitaxially grown on silicon wafers 
and allows one to engineer the bandgap, energy band structure, effective masses, mobilities 
and numerous other properties while fabricating circuits using conventional Si processing 
and tools.i^ 

There has been a tremendous effort in exploring the ways of breaking the silicon lattice 
symmetry and mixing different momentum states in order to induce optical gain. 3,4 The 
radiative efficiency depends on the competition between non-radiative fast processes and 
relatively slow radiative processes. To optimize the efficiency we have to eliminate non- 
radiative channels by using high purity materials and increase the oscillator strength of 
radiative channels. The above criteria also apply to photovoltaics. One modality for im- 
proving the optical response is quantum confinement. Confinement of the charge density in 
quantum wells (QWs) permits the relaxation of the optical selection rules for an interband 
transition. In addition, the band folding in superlattice (SL) structures will enhance the 
absorption. Consequently a Si based SL will reduce the symmetry which translates into 
band folding toward the zone center; and as a result vertical transitions will be available at 
energies closer to the indirect bandgap. The discovery of visible light emission from porous 
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FIG. 1: The band structure of bulk Si with the parameters of lVogl et al 
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adapted with the spin- 



orbit constant A = 0.045 eV. The single group notations are used since the spin-orbit splitting 
cannot be distinguished on the chosen scale. 



silicon has enticed many researchers to reactivate studies of the optical properties of silicon 
based nanostructures£ Efficient visible light emission has been observed in Si-CaF2 SLs, 6 
Si-SiC>2 SLsji Si nanopillar structures^ and in Si crystallites £ 

Beryllium chalcogenides have attracted much attention for the blue wavelength of region 
light emitting devices or UV detectors . 10 ^ 1 They provide long life-time II- VI lasers because 
of the strong covalent bonding of Be. BeTe and BeSe are wide-bandgap zinc blende semicon- 
ductors with lattice constants on either side of Si: they have the lattice constants of 5.6269 
and 5.1477 A, respectively, 3.6 % larger and 5.2 % smaller than Si. Vegard's law indicates 
that the lattice matched composition with Si is BeSeo.41Teo.59. Thus recent developments^*^ 
in the growth of silicon lattice- matched BeSe .4iTe .59 open the opportunity for a new class 
of Si based devices. 

In this paper we analyze quantum confinement and band folding in silicon based het- 
erostructures (Si/BeSeTe heterostructures) as ways of enhancing light absorption. A 
semiempirical tight-binding (TB) method^ is used throughout the study. Although nowa- 
days ab-initio methods (as they are formulated within density-functional theory)^ are be- 
coming routinely utilized for electronic structure and optical properties, they underestimate 
the bandgaps and therefore they need to be used with care. Moreover, the bandgap problem 
is solved within the so called GW ( G = Green function, W = screened Coulomb interaction) 



approximation which is a computational demanding task and rarely used for large systems 
like ours (in our case there is an alloy plus a heterostructure) .— In this work we choose a 
TB method together with a new method^ intended to deal with highly mismatched semi- 
conductors, thus considering band-bowing of the BeSeTe alloy This paper is organized as 
follows. In section II we discuss methods to improve the light absorption in Si based het- 
erostructures such as quantum confinement in a single (uncoupled) QW and band folding in 
SL structures. A significant increase in direct absorption is obtained for a narrow uncoupled 
silicon QW. In this spirit, we calculate the electronic and optical properties of Si/BeSeTe 
SLs using the TB method mentioned above. In the last section the conclusions are outlined. 



II. IMPROVED OPTICAL ABSORPTION BY QUANTUM CONFINEMENT 
AND BAND FOLDING 



A. Band Structure of BeSeo.4iTeo. 
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Since both Si and BeSe .4iTe .59 (as we will see below) are indirect band materials, full 
band calculations are required for a theoretical understanding of electron transport and 
optical properties in such heterostructures. For this purpose we use the empirical tight- 
binding (ETB) method as one of the most used tools in research of complex molecular and 
solid-state systems.^ Despite the fact that ETB is based on physical approximations such 
as the one-particle picture, short-range interactions, etc. (a posteriori justified, however), 
it gives fast and satisfactory results. The short-range nature of the model is suitable for 
modeling heterointerfaces which are present in such quantum structures. One of the most 
popular TB models is the sp 3 s* model.— This is a twenty-band model if spin-orbit coupling 
is included.— The presence of an additional s* orbital is able to reproduce th e bandgap i n 



indirect semiconductors like silicon. The TB para meters o f silicon are those from lVogl et al 
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augmented with spin-orbit coupling according to 

BeTe and BeSe are quite new materials in the sense that there are few experimental 
facts about these semiconductors. The bulk TB parameters of BeTe and BeSe were deter- 
mined by fitting the GW calculations of Fleszar and Hanked The band bowing has been 
included through a generalization of the band anti-crossing model for highly mismatched 



semiconductors.— The details of the model and the TB parameters can be found in Ref. 
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FIG. 2: Energy bands for BeSeo.41Teo.59 calculated with TB parameters from Ref. Il8t 

The energy band diagram for BeSeo.41Teo.59 is shown in Fig. |2l The material is seen to 
be indirect bandgap with the conduction band minimum at the X point. The predicted 
fundamental bandgap is 1.83 eV. The direct bandgap is 3.59 eV and the spin splitting is 
0.67 eV. Previous electrical measurements^^ indicate a conduction band offset of at most 
1.2 eV for the Si/BeSe .4iTe .59 heterostructure. As we will see in the following subsections 
there are interface subbands induced at the heterostructure interfaces. This fact may over- 
estimate the conduction band offset. Therefore we choose a value of 1 eV for the conduction 
band offset. This makes the Si/BeSeo.4iTeo.59 a type II heterostructure. Thus the top of the 
valence band in Si is higher than in the top of the valence band in BeSeo.41Teo.59. The TB 
parameters and the conduction band offset will be used in the following sections in order to 
calculate electronic and optical properties of the Si/BeSe .4iTe .59 SL. 

B. A Single Quantum Well 

The absorption edges in bulk silicon are indirect and the jumping of the electron from 
valence band (VB) to conduction band (CB) is assisted by a phonon (in electron transitions 
from VB to CB the k momentum must be conserved). The band structure of Si is shown 
schematically in the Fig. [TJ The indirect bandgap is Ai — r„ 5 J. The absorption process 
involves two steps, one is electron-photon interaction and the other is electron-phonon in- 
teraction. Basically we have two processes denoted by 1 and 2. In process 1 an electron is 
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first excited to the Tis state and then by an emission of a phonon the electron arrives in the 
Ai CB state. Similarly for process 2, however, in this process it starts with an emission of 
a phonon. For this reason, there is a long absorption tail between 1.12 eV and and about 3 
eV that reflects the indirect nature of the bandgap. Then a sharp rise in absorption occurs 
starting around 3.2 eV (380 nm) that is associated with the direct transition at T point 
(T 25 5 — > Ti 5 ) whose energy is 3.4 eV (365 nm). 

In direct bandgap materials light absorption in 2D systems is formally similar to 3D 
systems unlike the light absorption in indirect bandgap systems. The absorption coefficient 
for an uncoupled Si/ BeSeo.41Teo.59 QW has basically three components^ 



a M = A 



Pth$ (Eg + E e + E h -hw)+J2 Peh ^x q $ {E g + E e + E h - Huj - hw Xq ) 



eh,X 



•i 



+ + l ) 6 ( E 9 + E e + E h -hLJ + Hu Xq ) 



,(1) 



where A is a constant, pf h gives the direct bandgap contribution and p^ h gives the phonon 
assisted contribution with 1-phonon absorption, and p e eh is the phonon assisted contribution 
with 1-phonon emission. Basically, 

P d eh = Pcvhh (ko) , (2) 
Plh,\q = PehM = Va,Ieh (q - k ) R X . (3) 

p cv is the bulk dipole matrix element between bands, R\ is the matrix element contribution 
from electron-phonon interaction, and I e h is the overlap between the envelope-functions at 
different Brillouin points: 

/efc(q) = J dr^(r)Vh(r) e~^\ (4) 

where ko is the location in A;— space of each conduction valley. For silicon, if z is the growth 
direction, then the 2 valleys along z-directions (with k y = k x = 0) are responsible for direct 
transitions. The other four valleys contribute to phonon-assisted absorption. For infinite 
wells the overlap integral I e h has the behavior depicted in Fig. El For very narrow QW's the 
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FIG. 3: The overlap integral (Eq.Q) for infinite wells in silicon. The wave vector ko in Eq.(£Q) is 
the wave vector for silicon conduction valleys along z— axis. 

overlap integral tends to one, which means that in a genuine 2D system only the transverse 
momentum has to be conserved. For the other asymptotic limit, i.e. very wide well, the 
overlap integral vanishes. In physical terms, this says that the system became genuinely 3D 
and any electromagnetic transition has to be vertical in the absence of phonons. Moreover, 
from Fig. |3] we may expect to have a strong direct transition (with an overlap integral no 
less than 0.1) for a QW with a width up to 20 A, i.e. an ultrathin QW. In the following we 
will discuss the possibility of improving the oscillator strength in Si/ BeSeo.41Teo.59 based 
SLs. 

C. Electronic Structures of (Si2)n/(BeSeo.4iTeo.5g) m [001] Superlattices 

The use of SL structures relaxes the condition of ultrathin QW. The alternation of QWs 
and barriers along [001] growth axis will generate a band folding in the SL and band mixing of 
zone-center and zone edge states . 23 ' 24 The band mixing will inherently enhance the oscillator 
strength for direct transition in Si structures. On the other hand, band folding will induce 
states in which vertical transitions are possible at energies far lower than 3.4 eV, the lowest 
energy for vertical transitions in bulk Si. 

We consider a Si/BeSeo.4iTe .59 SL whose layers are perpendicular to [001] direction. We 
employ the nearest neighbor sp 3 s* TB Hamiltonian including spin-orbit interaction. Al- 
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(Si 2 V(BeSe M1 Te J 2 (Si 2 ) 8 /(BeSe 041 Te a59 ), 




though there are optimized nearest-neighbor TB parameters in the literature 2 ^ we chose the 
parameters given in Ref . [lj. These parameters include the chemical trends which are needed 
for an overall description of both valence and conduction bands. As it was noted before, 
the optimized nearest-neighbor parameters may have unphysical interaction energies and 
they are dedicated for either valence band or conduction band calculations.— Moreover the 
nearest-neighbor model has the inability to reproduce the transverse/light X-valley effective 
mass. However, the conduction band edges in low- dimensional structures will be determined 
by the longitudinal/heavy effective mass at the conduction band minimum. Therefore we 
opted for a less accurate yet a very representative set of parameters because these parameters 
are within the physical capabilities of the nearest neighbor model. 

We denote this SL as a (Si 2 )„/(BeSeo.4iTe .59) m SL with n two-atom thick layers of Si 
and m two- atom thick layers of BeSeo.41Teo.59 repeated periodically. The defined supercell 

consists of 2 (n + m) adjacently bonded atoms Si, Si, Si, Be, Te/Se, Be, Te/Se, Be, 

Te/Se. Covalent bonds between Si on the one hand and Be or Se/Te on the other are not 
electrically neutral.— For Si-based heterostructures special care has to be exercised due to 
excess/deficiency of electrons at interface and due to antiphase domains. 27 These problems 
can be overcome by vicinal growth and passivated surfaces . 12113 In this way one can minimize 
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the interface charging effects. Moreover, self-consistent ab-initio calculations show that the 
electric fields due to difference in ionicity between the two semiconductor constituents of 
a superlattice dominate the total electric field.— Therefore the accumulated charge with 
respect to the bulk contributes in a much lesser degree to the total electric field in the 
superlattice structure. In addition to that, the ETB fully self-consistent calculations as 
those described by Delia Sala and coworkers^ are not much different from those non-self- 
consistent calculations provided that the current densities in the confined system are low.— 
All the above facts lead us to conclude that a "zero- field" approximation will give us at least 
a qualitative picture of Si/BST structures. 

The TB matrix elements of the Si/BeSeo.4iTe .59 SL are taken directly from those bulk 
values. The on-site energies of the BeSe .4iTe .59 are accordingly changed to match the 
valence band offset at the interface. Simple averages were used to supply the parameters 
connecting different materials at the interface. Since the spin orbit must be included, the 
SL Hamiltonian will be represented as having 20 (n + m) functions. Once the TB matrix 
elements have been established the SL band structure reduces to the diagonalization of the 
20 (n + m) x20 (n + m) matrix Hamiltonian. The band structure of (Si2) n / (BeSeo.4iTeo.59) m 
SL for m = 2 and n = 6, m = n = 4, m = n = 8, and m = n = 6 are displayed in Fig. 01 The 
zero of energy corresponds to the top of the valence band in bulk Si. The band folding effect 
can be seen as many crowded subbands. Two interface bands (l v and I c ), one empty and one 
occupied, were found. They lie in the lower and upper parts of the bandgap of bulk silicon, 
respectively. The origin of these interface bands rests on the polar nature of the interface as 
was also found in GaAs/Ge SL . 31 i 32 i 33 The polarity of the interface originates from the large 
differences in the on-site energies for the constituent atoms (Si and Be or Se/Te). Even if a 
(110) non-polar interface is used, one interface band is still found in II-VI/IV SLs . 34 ' 35 We 
calculated the planar charge density of some of the occupied and empty band edge states 
of the SL with m = n = 4 at the T point in the Brillouin zone. The planar charge density 
of some of the band edge states and interface states are depicted in Fig. El By interface 
states we mean states that die away within a few layers from the interface. We denote by 
T Iv , the interface states of the I v band and by T Ic , the interface states of the I c band at V 
point. The charge density of the T Iv interface state has a maximum at the Be-Si interface, 
while the next occupied state is more confined in the BeSeo.41Teo.59 layer, but it still spills 
over into the Si barrier. The third occupied state spills more over into the Si slab, not only 
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FIG. 5: Planar charge densities in the (Si2)4/(BeSeo.4iTeo.5g)4 [001] superlattice at T point. The 
occupied states are the interface state T Iv and the first top two confined states (T V1 and T V2 ). 
The empty states are the interface state T Ic and the first two confined states (T C1 and T C2 ). The 
solid vertical lines denote the interfaces. 

due to interface but also due to the fact that the state is closer to the top of hole barrier in 
Si. The charge density of the T Ic interface state has a maximum at the Si-Se/Te interface, 
while the next empty state is confined, but it is still localized toward the interfaces. Only 
the third empty state is genuinely confined to the silicon slab. We believe that the charge 
distribution of the interface subbands is related to the acceptor behavior of Be and donor 
behavior of Se and Te with respect to Si. 36 



D. Optical Properties of (Si2) n /(BeSeo.4iTeo.59) m [001] Superlattices 

The electronic contribution to absorption spectrum is given by a a b s (u>) ~ u e 2 (to), where 
Cabs (w) is the absorption coefficient and^ 

2n 2 he 2 

£ 2 (oj) = — } f CV)k 5 (E c k - E v k - hu) (5) 

c,v,k 

is the imaginary part of the dielectric function. Here m is the electron mass, Q is the vol- 
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ume, e is the electron charge, H is the Planck constant and f cv ^ is the oscillator strength for 
the direct transition from the state |i>, k) to |c, k), with the photon momentum neglected. 
In Eq. (j5J) it is invoked the electric dipole approximation in the low-temperature case. In 
ETB the expressions for optical absorption and, in general, the interaction with the elec- 
tromagnetic field are different than those encountered in conventional quantum mechanics 
due to incompleteness of the basis.— In the electric dipole approximation, however, the ex- 
pressions of optical absorption are the same. The full temperature-dependent equations for 



optical absorption are derived in Refs. |38| and 39. Beyond the dipole approximation the 
standard expression (JHJ) is no longer valid. The problem is treated exhaustively in Ref. Q, 
where it is explicitly shown how the Hamiltonian is altered by the vector potential of the 
electromagnetic field. 

The oscillator strength is defined as 

_ 2 |(c,k| e-p \v,k}\ 2 

Jcv,k 7—1 7-i Y°) 

m &c,k — &v,k 

In Eq. (JHJ) \v, k) and |c, k) are the valence and conduction band eigenstates, E v ^ and E c ^ are 
their corresponding energies, e is the polarization of light, and p is the momentum operator. 
In the empirical tight-binding approach the momentum matrix element is defined as^ 

Tfl 

(c, k| p \c, k) = - (c, k| V k tf (k) \v, k) . (7) 

h 

Although Eq. ((Zj) does not consider the intra-atomic part of the momentum matrix 
element,— i 4 - 2 - it is gauge invariant . 38 ' 39 Its accuracy can be improved by going beyond near- 
est neighbor interaction or by using non-orthogonal orbitals.— We first calculate the joint 
densities of states (JDOS) because the variation of the oscillator strength over the Brillouin 
zone is small.— We assume that light is propagating along the SL growth direction. The 
JDOS represents the number of states that can undergo energy and k-conserving transitions 
for photon frequencies between u and uj + du. The JDOS associated with Eq. Q are shown 
in Fig. IHlfor SLs with n = 12 and m = 4, n = m = 8, n = 8 and m = 4, and m = n = 6. A 
0.05 eV broadening was considered for each electronic energy. The summation over Brillouin 
zone was replaced by the summation over special points in the Brillouin zone . 44145 

Due to band folding the absorption edges for vertical transitions are lowered toward the 
indirect bandgap of bulk Si. Moreover, the absorption edges are not determined by the 
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interface subbands but by the energy difference between the lowest empty confined state 
and highest occupied confined state. Thus the absorption edges for vertical transitions are 
below the indirect bandgap of Si because the heterostructure is of type II. The curves rise 
slowly and, with increasing m(n), the absorption edges extend to lower energies. However, 
the thickness of the BeSe .4iTe .59 layers play a greater role than the thickness of Si layers 
in the variation of the absorption edge. It can be explained by the difference of the electron 
effective masses in Si and hole effective masses in BeSeo.41Teo.59, i.e., the heavy hole effective 
mass in BeSeo.41Teo.59 is smaller than the effective mass of heavy electron valley in Si. Simi- 
lar results were found in the calculations for porous Si with periodic boundary conditions?^ 
the absorption band edge does not vary much with Si porosity. The strength of the optical 
absorption is also determined by the oscillator strength, thus we can also check if the tran- 
sition are allowed or not allowed. We denote by Iv, VI, and V2 the first three top valence 
subbands, and by Ic, CI, and C2 the first three conduction subbands. We calculated the 
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oscillator strengths of several interband transitions relative to the oscillator strength of the 
direct transition in bulk Si. The results are shown in Fig. [7|for m = n = 4. The oscillator 
strengths of interband transitions are at least 10 times smaller than their bulk counterpart 
and basically range from 10 ~ 3 to 10 _1 relative to r 25 > — > T 15 transition in bulk Si. The 
strongest transitions are those coming from the interface subband and those coming from 
the first confined hole level to the confined electron levels. 

III. CONCLUSIONS 

We studied the electronic and optical properties of silicon-based heterostructures 
(Si/BeSeo.4iTeo.5g). Based on the fact that significant light absorption occurs for thin uncou- 
pled Si quantum wells we studied the (Si) n /(BeSeo.4iTeo.5g) m superlattice structures with a 
semiempirical tight-binding method in the sp 3 s* nearest neighbor model including spin-orbit 
interaction. The tight-binding model considered large band bowing in the BeSeo.41Teo.59 al- 
loy. The heterostructure Si/BeSeo.4iTeo,59 was found to be of type II. Electronic structure 
and optical properties were calculated for various superlattice structures. Two interface 
bands were found to exist in the bandgap of bulk silicon. We found that band folding 
induces vertical transitions below the indirect bandgap of bulk Si. In addition, calculated 
oscillator strengths for vertical transitions near the optical band edge show the mixing of the 
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zone-center and zone edge states of bulk Si for conduction subbands. Therefore, the transi- 
tions are optically allowed and the response of silicon based heterostructures to illumination 
is enhanced. 
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